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Abstract 

Organophosphorus  (OP)  nerve  agents  are  potent  suicide  inhibitors  of  the  essential  neurotransmitter-regulating  enzyme 
acetylcholinesterase.  Due  to  their  acute  toxicity,  there  is  significant  interest  in  developing  effective  countermeasures  to  OP 
poisoning.  Here  we  impart  nerve  agent  hydrolysis  activity  into  the  human  drug  metabolism  enzyme  carboxylesterase  1. 
Using  crystal  structures  of  the  target  enzyme  in  complex  with  nerve  agent  as  a  guide,  a  pair  of  histidine  and  glutamic  acid 
residues  were  designed  proximal  to  the  enzyme's  native  catalytic  triad.  The  resultant  variant  protein  demonstrated 
significantly  increased  rates  of  reactivation  following  exposure  to  sarin,  soman,  and  cyclosarin.  Importantly,  the  addition  of 
these  residues  did  not  alter  the  high  affinity  binding  of  nerve  agents  to  this  protein.  Thus,  using  two  amino  acid 
substitutions,  a  novel  enzyme  was  created  that  efficiently  converted  a  group  of  hemisubstrates,  compounds  that  can  start 
but  not  complete  a  reaction  cycle,  into  bona  fide  substrates.  Such  approaches  may  lead  to  novel  countermeasures  for  nerve 
agent  poisoning. 
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Introduction 

Organophosphorus  (OP)  nerve  agents  are  some  of  the  most 
poisonous  chemicals  known  (Figure  1A)  [1],  By  covalently 
phosphonylating  the  catalytic  serine  of  the  neurotransmitter¬ 
regulating  enzyme  acetylcholinesterase  (AChE),  OPs  produce  a 
cholinergic  crisis  that  can  cause  respiratory  failure  and  death 
[2].  OP  nerve  agents  have  been  used  as  chemical  weapons 
by  both  established  governments  and  terrorist  groups  [2,3]. 
In  addition,  there  are  an  estimated  250,000  deaths  world¬ 
wide  each  year  associated  with  OP  pesticides  [4],  Thus,  there 
is  significant  interest  in  developing  novel  approaches  to 
detoxify  these  compounds  for  military,  security  and  clinical 
applications  [5], 

Human  carboxylesterase  1  (hCEl)  is  primarily  expressed  in 
the  liver  and  metabolizes  a  variety  of  chemically  distinct 
substrates,  from  potentially  harmful  endogenous  compounds  to 
exogenous  compounds  including  drugs  and  environmental  toxins 
[6].  hCEl  is  homologous  in  structure  and  catalytic  mechanism 
to  AChE  and  exhibits  features  that  make  it  attractive  as  a 
potential  nerve  agent  hydrolase  [6].  Both  hCEl  and  AChE 
employ  a  catalytic  triad  consisting  of  serine,  histidine,  and 
glutamic  acid  residues.  In  addition,  for  both  enzymes,  OP  nerve 


agents  act  as  hemisubstrates  -  compounds  that  can  start  but  not 
complete  a  reaction  cycle  to  regenerate  free  active  enzyme 
[7,8,9].  Indeed,  hCEl  is  potently  inhibited  by  nerve  agents, 
forming  long-lived  covalent  phosphonyl-enzyme  complexes 
(Figure  IB)  [8].  We  have  previously  shown  that  hCEl  does 
not  undergo  the  dead-end  “aging”  process  from  this  interme¬ 
diate  that  is  observed  for  the  cholinesterases,  AChE  or 
butyrylcholinesterase  (BChE)  [10].  In  addition,  wild-type  (wt) 
hCEl  exhibits  a  slow  hydrolysis-based  reactivation  after 
exposure  to  sarin,  but  no  reactivation  is  observed  when  hCEl 
is  exposed  to  soman  or  cyclosarin  [8]. 

Thus,  we  sought  to  design  an  hCE  1  variant  capable  of 
hydrolyzing  OP  nerve  agents.  We  employed  molecular  modeling 
techniques  based  on  crystal  structures  of  hCEl  in  trapped  covalent 
complexes  with  agents  to  engineer  amino  acids  capable  of 
promoting  dephosphonylation.  Herein,  we  report  the  placement 
of  a  pair  of  residues  in  the  hCEl  active  site  that  significandy 
increases  the  rate  of  nerve  agent  hydrolysis  relative  to  native 
hCEl,  while  still  retaining  high  affinity  and  selective  functional 
substrate  activity  towards  these  compounds.  Thus,  using  only  two 
amino  acid  changes,  we  have  created  one  of  the  most  efficient 
enzymes  engineered  to  date  with  respect  to  converting  hemisub¬ 
strates  into  substrates. 
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Figure  1.  Organophosphate  (OP)  inhibition  of  human  carboxylesterase  1  (hCEl).  A.  Three  G-type  OP  nerve  agents  and  OP  model 
compound  (R  represents  respective  O-alkoxy  groups).  Wild-type  hCEl  preferentially  binds  the  stereoisomers  shown  (7).  B.  Schematic  mechanism  of 
OP  hydrolysis  by  hCEl.  X  represents  the  leaving  group  and  *  denotes  a  non-reactive  state. 
doi:1 0.1 371  /journal.pone.001 7441  .g001 


Results 

Structure-Guided  Dyad  Design 

To  hydrolyze  the  covalent  phosphonyl-enzyme  OP-bound 
species  of  hCEl,  we  introduced  a  histidine  (H)  at  position  146 
and  a  glutamic  acid  (E)  at  position  363.  These  residues  were 
chosen  because  they  appeared  ideally  positioned  to  order  a  water 
molecule  for  hydrolysis.  The  corresponding  wild-type  residues 
VI 46  and  L363  are  positioned  on  either  side  of  the  active  site 
pocket,  adjacent  to  the  OP-serine  bond  formed  during  the  two- 
step  catalytic  mechanism  of  the  enzyme  (Figure  2).  VI 46  is  part 
of  an  extended  loop  between  (37  and  a2  that  also  contains  the 
oxyanion  hole  formed  by  residues  G142  and  G143.  L363  caps 
a  12,  at  the  edge  of  the  hCEl  active  site.  As  outlined  below,  the 
structurally  similar  residue  glutamine  (Q)  was  used  as  a  control 
mutation  for  each  position.  hCEl  mutants  were  successfully 
expressed,  purified,  and  shown  to  exhibit  enzyme  activity  using  a 
standard  esterase  substrate  para-nitrophenylbutyrate  (pNPB)  [1 1], 
although  wildtype  hCEl  was  ~30-fold  faster  toward  pNPB 
than  the  V146H/L363E  mutant  (660  nmol/min/mg  vs.  21 
nmol/ min/ mg). 

Enhanced  Nerve  Agent  Hydrolysis 

The  V146H/L363E  hCEl  mutant  exhibited  a  significant 
increase  in  enzyme  reactivation  following  nerve  agent  inhibition 
(Figure  3A).  We  had  shown  previously  that  wild-type  (wt)  hCEl 
reactivates  following  sarin  exposure  at  a  rate  (A3)  of 
2.4±0.50xl0  4  min  ,  or  a  half-time  of  reactivation  (ti/2)  of 
46  hrs,  and  that  wt  enzyme  was  not  able  to  recover  following 
soman  or  cyclosarin  inhibition  [8].  In  contrast,  however,  the 
V146H/L363E  hCEl  mutant  spontaneously  reactivated  following 
sarin  and  soman  inhibition  with  rates  of  12±2.0xl0-  min- 
(ti/2  °f  9.5  hrs)  and  10±  1.0  x  10-4  min-1  (t1/2  of  11.5  hrs), 


respectively.  The  reactivation  rates  of  V146H/L363E  hCEl  after 
sarin  and  soman  inhibition  were  5-  and  20-fold  faster,  respectively, 
than  the  wt  enzyme.  Most  strikingly,  the  designed  V146H/L363E 
dyad  exhibited  a  rate  of  reactivation  of  67±6.0  x  10  4  min  1  (ti/2 
of  1 .2  hr)  following  cyclosarin  exposure,  an  increase  in  reactivation 
of  ~  33,000-fold  relative  to  wt  enzyme  (where  the  wt  rate  was 
estimated  to  be  1  x  1 0  5  min  1 ,  the  lower  limit  of  rate  measure¬ 
ment).  Recall  that  no  reactivation  was  observed  with  wt  hCEl  in 
the  presence  of  this  OP.  Thus,  efficient  recovery  of  carboxylester¬ 
ase  activity  following  nerve  agent  exposure  was  achieved  by 
introducing  a  pair  of  residues  adjacent  to  the  hCEl  active  site. 

We  next  tested  a  series  of  single  and  double  mutant  controls 
containing  Q_  residues  at  positions  VI 46  and  L363  to  examine  the 
effect  of  each  mutant  on  reactivation  rate  (A3)  with  respect  to  sarin, 
cyclosarin,  and  soman  hydrolysis  (Figure  3B,  Table  1).  For  sarin 
hydrolysis,  introduction  of  a  polar  side  chain  (H  or  Q)  in  either 
position  alone  or  in  combination  produced  similar  rates  of 
reactivation,  ranging  from  3.2  xlO-4  to  26xl0-  min-  ,  up  to 
a  11 -fold  increase  relative  to  wt  hCEl.  Following  soman  or 
cyclosarin  inhibition,  L363E  appeared  critical  to  increasing  A3 
from  sub-detection  limits  in  wt  hCEl  to  4.0  xlO-  and 
19x10  4  min  ',  up  8-  and  9,300-fold,  respectively  (Table  1). 
For  soman,  polar  residues  at  position  146  (H  or  Q)  weakly 
stimulated  reactivation;  however,  when  V146FI  was  combined 
with  L363E,  reactivation  increased  20-fold  relative  to  wt  hCEl 
(Table  1). 

However,  the  largest  A3  enhancement  was  for  cyclosarin  and  the 
V146H/L363E  mutant.  This  designed  enzyme  exhibited  an 
increased  rate  of  enzyme-mediated  cyclosarin  dephosphonylation 
of  >33, 000-fold  relative  to  the  wt  enzyme  (Table  1).  A  second 
double-mutant  in  which  a  Q^at  position  146  was  combined  with  an 
E  at  position  363  was  also  found  to  exhibit  robust  enzyme 
reactivation  following  cyclosarin  exposure,  with  an  rate  increase  of 
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Figure  2.  Human  carboxylesterase  1  active  site  structure.  Active  site  of  human  carboxylesterase  1  covalently  inhibited  via  S221  with  cyclosarin 
(magenta)  [8],  The  other  catalytic  residues,  in  addition  to  S221,  are  H468  and  E354  (yellow),  and  are  surrounded  by  hydrophobic  residues  (grey 
surface)  including  VI 46  and  L363  (light  blue),  as  well  as  the  oxyanion  hole  (white). 
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>21, 000-fold  relative  to  wt  hCEl  (Table  1).  Thus,  a  glutamic 
acid  at  position  363  appeared  to  function  as  a  catalytic  base  to 
enhance  cyclosarin  reactivation;  a  glutamine  in  that  position  failed 
to  stimulate  enzyme  activity  in  the  presence  of  cyclosarin.  Taken 
together,  these  data  demonstrate  that  targeted  mutagenesis  within 
hCEl  can  significantly  increase  nerve  agent  processing  by  this 
serine  hydrolase. 

Structural  Mechanism  for  Enhanced  Cyclosarin  Hydrolysis 

The  crystal  structure  of  hCEl  in  covalent  complex  with 
cyclosarin  was  then  used  to  generate  an  energy-minimized  model 
of  V146H/L363E  hCEl.  In  the  resultant  model,  the  side  chains  of 
L363E  and  V146EI  were  found  to  be  3.3  A  apart.  In  addition,  the 
orientation  of  L363E  suggested  that  a  water  molecule  could  be 
positioned  3.2  A  from  L363E  and  3.2  A  from  the  phosphorous 
atom  on  the  OP-enzyme  adduct  (Figure  4A).  Based  on  this  model 
and  the  biochemical  data  presented  above,  we  propose  that  L363E 
is  positioned  for  water  activation  and  is  stabilized  by  HI 46 
(Figure  4B).  Indeed,  pH  rate  dependence  studies  showed  that 
reactivation  proceeded  most  effectively  below  pH  6.2,  when 
V146H  is  likely  to  be  positively  charged  (Figure  4C). 

Bimolecular  Rates  of  Inhibition 

To  ensure  that  the  enhanced  rates  of  reactivation  of  V146H/ 
L363E  were  not  artifacts  of  altered  substrate  binding,  we 
determined  the  bimolecular  rates  of  inhibition  (£;)  of  the  mutant 
protein  with  racemic  cyclosarin  or  stereoisomers  of  OP  model 
compounds  similar  in  structure  to  cyclosarin,  and  compared  those 
rates  to  wt  enzyme.  V146H/L363E  hCEl  exhibited  a  similar  A-  to 
wt  hCE  1 ,  confirming  that  the  ability  of  the  mutant  protein  to  bind 


cyclosarin  and  related  OP  compounds  was  not  disrupted  (Table  2, 
Table  3).  Indeed,  as  measured  with  stereoisomers  of  the  model 
compounds,  V146H/L363E  exhibited  an  improved  binding 
constant  (A^)  compared  to  wt  hCEl  with  either  stereospecific 
compound  [8].  Therefore,  the  addition  of  an  ion  pair  across  the 
active  site  entrance  was  not  detrimental  to  cyclosarin  binding. 

Discussion 

OP  nerve  agents  serve  as  hemisubstrates  of  the  wild-type  AChE, 
BChE,  and  hCEl  enzymes  because  these  chemicals  fail  to 
complete  the  typical  two-step  serine  hydrolase  mechanism. 
Rather,  OP  nerve  agents  produce  relatively  long-lived  covalently 
phosphonylated  active  site  serines  [12,13].  Histidine  (H468  in  wt 
hCEl;  see  Figure  2B)  activates  the  serine  (S221  in  hCEl)  for 
nucleophilic  attack  on  the  OP  phosphonyl  center.  The  resultant 
pentahedral  intermediate  is  stabilized  by  hCEl  backbone  nitrogen 
atoms  in  the  oxyanion  hole  (G142  and  G143);  the  fluoride  leaving 
group  dissociates  upon  collapse  of  the  transition  state,  resulting  in 
a  covalently  modified  enzyme  [14].  The  resulting  tetrahedral 
phosphonyl  adduct  adversely  affects  the  ability  of  H468  to 
facilitate  base  mediated  C-O  bond  cleavage  [15],  which  leaves 
this  residue  in  a  state  where  it  cannot  activate  the  water  molecule 
required  for  hydrolytic  desphosphonylation  [16].  Thus,  for  serine 
hydrolases  like  AChE  and  hCEl  that  are  inhibited  by  OP 
hemisubstrates,  the  addition  of  a  strong  oxime  is  often  required  to 
complete  the  catalytic  cycle  and  regenerate  active  enzyme. 
Current  treatments  following  OP  exposure  include  administration 
of  strongly-nucleophilic  oximes  to  dephosphonylate  AChE. 
However,  such  compounds,  such  as  2-pralidoxime  obidoxime, 


PLoS  ONE  |  www.plosone.org 


3 


March  2011  |  Volume  6  |  Issue  3  |  el  7441 


Enzyme  Design  for  Novel  Function 


A) 


120 


100 


~  80 


.=  60 


c 

'c 

D 


Cyclosarin 


i-— 1 


Sann 


* 

TJ 


VI46H/L363E  hCEl 


/# 

/  * 

.-i 


Soman 


Sarin 


40 


20 


y 

I 


Wild  type  hCEl 


i 


Soman 


Cvclosarin 


0  10  20  30  40  50  60  70  80 

Time  (hrs.) 


■Sarin 
®  Soman 
■  Cyclosarin 


Figure  3.  Reactivation  of  hCEl  following  nerve  agent  exposure.  A.  Spontaneous  reactivation  of  V146H/L363E  hCEl  following  inhibition  by 
racemic  sarin  (blue),  soman  (green),  or  cyclosarin  (red).  Wild  type  hCEl  (grey)  only  reactivates  following  sarin  inhibition  (7).  n  =  6,  s.d.  B.  Rates  of 
dephosphonylation  for  hCEl  variants  in  the  presence  of  sarin  (blue),  soman  (green)  and  cyclosarin  (red).  n  =  3,  s.d. 
doi:1 0.1 371/journal.pone.OOI  7441  .g003 


do  not  offer  broad-spectrum  protection  against  all  agents  and  must 
be  administered  quickly  [17]. 

Based  on  previous  crystal  structures  of  hCEl  in  covalent 
complexes  with  nerve  agents  [7,8],  we  hypothesized  that 
introducing  a  general  base  catalyst  would  facilitate  the  activation 
of  a  properly  located  water  molecule  to  hydrolyze  an  OP  adduct. 
The  V146H/L363E  amino  acid  substitutions  dramatically  in¬ 
creased  the  spontaneous  rates  of  enzyme  reactivation  following 


sarin,  soman,  and  cyclosarin  inhibition  compared  to  native 
enzyme  (Figure  3).  A  panel  of  single  and  double  mutant  controls 
indicated  that  V146H  or  L363Q  facilitated  increased  sarin 
hydrolysis,  while  L363E  acted  as  general  base  catalyst  for 
cyclosarin  and  soman  hydrolysis.  Because  V146H,  L363Q,  or 
their  combination,  resulted  in  similar  rates  of  reactivation,  these 
mutations  may  enhance  dephosphonylation  by  interacting  with 
the  sarin  O-isopropyl  group,  thereby  affording  deprotonation  of 
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Table  1.  Rates  of  dephosphonylation  for  hCEl  mutants. 


hCEl 

Sarin 

Soman 

Cyclosarin 

k3  (xIO  4  min-1) 

Fold  Increase 

k3  (xIO  4  min-1) 

Fold  Increase  k3  (xIO  4  min  ’) 

Fold  Increase* 

WT 

2.4±0.5 

- 

0.5±0.3 

- 

<0.01 

- 

V146H 

20±4 

8 

2±1 

4 

0.032±0.006 

16 

V146Q 

3.2±0.9 

1 

2±1 

4 

0.34±0.09 

171 

L363E 

4.6  ±0.5 

2 

4±1 

8 

19±5 

9,300 

L363Q 

1 4±  1 

6 

<0.01 

- 

<0.01 

- 

V146Q/L363E 

5.7±0.7 

2 

2.6  ±0.5 

5 

43±7 

21,500 

V146H/L363Q 

26±3 

11 

<0.01 

- 

<0.01 

- 

V146Q/L363Q 

1 2±  1 

5 

<0.01 

<0.01 

V146H/L363E 

12±2 

5 

10±1 

20 

67±6 

33,355 

N.R.  is  no  reactivation  with  a  detection  limit  <1  xIO  5  min 
pH  7.4,  25 DC,  N  =  3,  s.d. 

*k3~1  xIO-7  min-1,  i.e.  (1  xIO-5  min-1)*(1%  activity)/100. 
doi:1 0.1 371/journal. pone.001 7441  .tOOl 


Figure  4.  Mechanism  of  reactivation  by  V146H/L363E  hCEl  after  cyclosarin  binding.  A.  Model  of  V146H/L363E  (cyan)  hCEl  with  PR 
cyclosarin  (magenta)  including  a  water  molecule  (red)  between  E363  and  the  central  phosphorus.  B.  Proposed  mechanism  for  enhanced  reactivation 
following  cyclosarin  inhibition.  C.  pH  dependence  of  V146H/L363E  (black)  and  L363E  (grey)  hCEl  dephosphonylation  following  cyclosarin  inhibition. 
doi:1 0.1 371/journal.pone.OOI  7441  .g004 
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Table  2.  Bimolecular  rates  of  inhibition,  Michaelis-Menten 
constants,  and  rates  of  reactivation  for  wild-type  and  V146H/ 
L363E  hCEl  against  racemic  cyclosarin  and  stereoisomers  of 
cyclosarin  model  compounds. 


V146H/L363E 

h  (xio3  M  ’s1 

Km  (xIO  6  M) 

k3  (xIO  4  min  n) 

Racemic3 

79±2 

N.D. 

67±6 

pBb 

5.4  ±0.4 

0.095  ±0.006 

50±10 

Psb 

0.12±0.02 

5.2±0.7 

40±10 

WT 

Racemic3 

11  ±5 

N.D. 

<0.0  T 

p»b 

1 5  ±  1 c 

0.32±0.01c 

<0.0  T 

Psb 

0.0053  ±0.0007c 

1 700±200c 

<0.0  T 

N  =  3,  s.d.,  N.D.  is  not  determined,  N.R.  is  no  reactivation,  pH  7.4,  25°C. 
aRacemic  bona  fide  cyclosarin, 
bstereoisomers  of  cyclosarin  model  compounds, 
c(8). 

doi:1 0.1 371  /journal,  pone.001 7441  .t002 

the  catalytic  H468  for  water  activation  and  hydrolysis.  A  G117H 
mutation  in  BChE  resulted  in  increased  rates  of  dephosphonyla- 
tion  through  an  analogous  mechanism  [16,18]. 

In  contrast,  the  L363E  mutation  positions  an  anionic  carbox- 
ylate  adjacent  to  the  phosphonyl-enzyme  intermediate  to  activate 
a  water  molecule  for  nucleophilic  attack.  Glutamic  acid  and 
aspartic  acid  residues  act  as  general  acid-base  catalysts  in  several 
established  enzyme  mechanisms,  including  those  of  lysozyme  and 
protein  tyrosine  phosphatases  [19,20],  and  apparently  support  the 
same  role  in  hCE  1 .  The  V 1 46H  addition  synergistically  increased 
base  activity  of  L363E  (Figure  4A).  The  pH  profile  for  cyclosarin 
hydrolysis  by  V146H/L363E  showed  that  the  fastest  rates  of  hCEl 
reactivation  was  when  HI 46  was  likely  to  be  in  a  charged  state, 
below  pH  6.2  and  where  a  cationic  HI 46  can  stabilize  the  anionic 
base  (Figure  4C).  Based  on  these  data,  we  postulate  that  E363 
Osl  is  stabilized  by  HI 46  and  Oe2  on  E363  is  available  to 
deprotonate  a  water  molecule  to  hydrolyze  the  covalent  cyclohexyl 
phosphonyl  group  (Figure  4B).  This  hypothesis  is  supported  by 
the  observation  that  the  V146Q/L363E  variant  also  exhibited  a 
significant  increase  in  cyclosarin  hydrolysis  (Table  1). 

Previous  attempts  to  introduce  general  base  catalysts  in  OP- 
inhibited  serine  hydrolases  have  resulted  in  enzymes  that  also 
showed  enhanced  rates  of  reactivations  against  soman,  sarin,  and 
VX,  but  the  ability  of  these  mutant  enzymes  to  bind  the  respective 
OPs  was  greatly  diminished  [10,16].  For  example  a  G117H/ 


E197(()  BChE  mutant  enhanced  the  rate  of  soman  hydrolysis  up  to 
2,500-fold  faster  than  wt  BChE,  but  the  soman  bimolecular  rate  of 
inhibition  for  this  mutant  was  decreased  by  a  similar  magnitude 
[10].  In  contrast  to  the  G1 17H/E197Q, BChE  mutant,  V146H/ 
L363E  hCE  1  exhibits  a  greater  second-order  rate  of  inhibition  (Aj) 
than  the  wt  hCEl  enzyme  (Table  2).  Further,  the  Km  values  for  the 
cyclosarin  model  compounds  reported  in  Table  2  demonstrate  that 
the  affinity  for  V146H/L363E  hCEl  is  greater  than  the  wt  enzyme 
for  these  analogs.  Taken  together,  these  data  show  that  V146H/ 
L363E  hCE  1  retains  greater  than  wt  affinity  for  cyclosarin,  and  that 
the  engineered  hCEl  mutants  enhance  rates  of  reactivation  via 
dephosphonylation  rather  than  decreasing  inhibition. 

In  conclusion,  we  have  rationally  designed  a  variant  form  of  the 
liver  detoxifying  enzyme  hCEl  that  spontaneously  dephosphony- 
lates  after  inhibition  by  sarin,  soman,  and  cyclosarin  up  to  33,000- 
fold  faster  than  wt  enzyme.  Wild  type  hCEl  has  a  catalytic 
efficiency  (£cat/7Tm)  [21]  towards  the  standard  esterase  substrate 
para-nitrophenylbutyrate  (pNPB)  of  1.2xl05  M  ’s  1  [22]. 

AChE,  one  of  the  most  efficient  enzymes  known,  exhibits  a  kCM/ 
Km  close  to  the  diffusion  control  limit  (108  M-1s-  )  for 
acetylcholine  hydrolysis  [23].  Combining  the  nanomolar  Am  of 
V146H/L363E  towards  the  P^  cyclosarin-like  compound  along 
with  rate  of  reactivation  after  inhibition  by  this  compound,  this 
double  mutant  shows  a  catalytic  efficiency  of  8.8x10'  M  s  *. 
Thus,  the  redesigned  hCEl  compares  favorably  to  other 
mammalian  enzymes  that  have  been  rationally  engineered  to 
improve  hemi-substrate  metabolism  (Table  4).  This  enzyme  will 
likely  require  substantial  increases  in  catalytic  efficiency  for  OP 
compounds  in  order  to  provide  in  vivo  protection,  but  nonetheless 
can  serve  as  a  lead  candidate  for  further  development  of  novel 
countermeasures  for  nerve  agent  or  pesticide  poisoning. 

Materials  and  Methods 

hCEl  Mutagenesis,  Expression,  and  Purification 

Site-directed  mutagenesis  was  done  using  PCR  to  introduce  the 
desired  amino  acid  changes  into  hCEl  (Genbank  accession 
M73499).  Briefly,  100  pM  of  both  sense  and  anti-sense  primers 
were  mixed  with  50  ng  hCEl  cDNA  in  a  pUC9  vector,  200  pM 
dNTPs,  1  x  Pfu- BSA  buffer,  and  2  u  Pfu  (New  England  Biolabs, 
Ipswich,  MA).  Mutations  were  produced  through  15  cycles  of 
PCR,  where  each  cycle  consisted  of  95°C  for  1  minute,  58°C  for 
30  seconds,  and  70  C  for  10  minutes.  The  reaction  was  digested 
for  1  hour  with  Dpnl  (Fermentas,  Burlington,  Ontario)  at  37°C 
and  transformed  into  chemically  competent  DH5a  cells  (Invitro- 
gen,  Carlsbad,  California).  pUC9  plasmids  were  isolated  with  a 
GeneJET  Plasmid  Miniprep  kit  (Fermentas)  and  mutations  were 
confirmed  by  DNA  sequencing.  Once  the  mutations  were 


Table  3.  Inhibition  and  Michaelis-Menten  constants  for  wild-type  and  V146H/L363E  hCEl  against  stereoisomers  of  sarin  and 
soman  model  compounds. 


hCEl 

Model  Compound 

k\  (xio3  m  V1) 

Km  (xIO  6  M) 

p« 

P  5 

P* 

Ps 

WT3 

Sarin 

0.77±0.05 

4.1  ±0.1 

1.4±0.2 

0.37±0.02 

Soman 

2.5±0.5 

0.001 5  ±0.0007 

0.22  ±0.02 

800  ±200 

V146H/L363E 

Sarin 

3.40  ±0.07 

4.7±0.2 

3.5±0.8 

o 

VO 

1  + 
O 

Soman 

10.8±0.4 

0.67±0.06 

0.101  ±0.005 

8±2 

(N  =  3,  s.d.) 
a(7). 

doi:1 0.1 371/journal. pone.001 7441  .t003 
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Table  4.  Catalytic  efficiencies  {kcat/Km)  of  engineered  enzymes  towards  hemisubstrates. 


Enzyme  Engineered  Activity  Mutation(s)  Catalytic  Efficiency  (M  1  s  ') 


Butyrylcholinesterase 

Sarin  hydrolysis3 

G117H 

0.79 

Glutathione  Transferase  Al -1 

Thioester  hydrolysisb 

A216H 

2.60 

4-Oxalocrotonate  Tautomerase 

Oxaloacetate  decarboxylase0 

P1A 

114 

Phospho(serine/threonine/tyrosine) 
Binding  Protein 

Phosphatase  activity0 

G120C 

490 

Human  Carboxylesterase  1 

Cyclosarin  analog  hydrolysis 

V146H/L363E 

877 

Cyclophilin 

Protease  activity0 

A9 1 S/F 1 04H/N 1 06D 

1675 

a(16). 

b(32). 

c(32). 

doi:1 0.1 371  /journal,  pone.001 7441  .t004 


successfully  incorporated,  the  1.7  kB  hCEl  cDNA  was  ligated  into 
pCIneo  for  mammalian  cell  expression  (Promega,  Madison,  WI). 
Non-secreted  forms  of  wt  and  mutant  hCEl  proteins  were 
expressed  in  COS-7  cells  (American  Type  Culture  Collection, 
Rockville,  Maryland)  as  previously  described  [11].  Presence  of 
active  hCEl  was  determined  by  measuring  the  rate  of  nitrophenol 
production  in  the  presence  of  3  mM  o-nitrophenol  acetate  through 
a  change  in  absorbance  at  420  nm.  Specific  activities  were 
determined  by  normalizing  hCEl  activity  relative  to  protein 
expression  as  visualized  through  western  blot  analysis  with  an  anti- 
hCEl  antibody  [11].  To  facilitate  measurement  of  additional 
kinetic  constants,  secreted  forms  of  wt  and  V146H/L363E  hCEl 
were  expressed  in  Spodoptera  frugiperda  Sf21  insect  cells  (Clontech, 
Palo  Alto,  California)  using  baculovirus-expression  vectors  and 
purified  as  previously  described  [24], 

Examination  of  Kinetics  of  hCEl  Variants 

Experiments  with  all  racemic  OP  nerve  agents  were  conducted 
at  the  United  States  Army  Medical  Research  Institute  of  Chemical 
Defense  (USAMRICD),  Aberdeen  Proving  Ground,  MD.  Nerve 
agents  were  obtained  from  the  U.S.  Army  Edgewood  Chemical 
Biological  Center  (ECBC,  Aberdeen  Proving  Ground,  MD). 
Analysis  by  nuclear  magnetic  resonance  spectroscopy  showed 
them  to  be  >95%  pure.  Dilute  nerve  agent  was  handled  according 
to  safety  guidelines  in  place  at  USAMRICD.  All  kinetic  assays 
were  conducted  at  25°C.  In  separate  experiments,  50  pL  of  whole 
cell  COS  lysates  expressing  an  hCEl  mutant  was  inhibited  with  a 
~  1 000-fold  molar  excess  of  racemic  sarin,  soman,  or  cyclosarin  for 
10  minutes.  Unbound  agent  was  removed  by  passing  inhibited 
samples  over  a  CENTRI-SEP  size  exclusion  column  (Princeton 
Separations,  Adelphia,  NJ).  The  column  eluate  was  diluted  10-fold 
in  0. 1  M  potassium  phosphate  buffer,  pH  7.4.  Aliquots  were  tested 
for  hCEl  carboxylesterase  functional  activity  using  5  mM  pNPB 
and  compared  to  an  uninhibited  sample.  Measurements  were 
taken  over  60  hours,  and  the  observed  rate  of  enzyme  reactiviation 
(Hobs),  or  the  ability  of  enzyme  to  hydrolyze  pNPB,  and  maximal 
percent  recovery  (Amax)  were  determined  by  fitting  the  collected 
data  to  equation  1: 

A=A0+Amax(l-e-kobs<)  (1) 

where  A  was  the  percent  activity  at  time  t,  and  A0  was  the  initial 
activity  at  #=0.  Reactivation  experiments  were  conducted  in 
triplicate  and  independently  replicated.  It  should  be  noted  that 
these  experiments  gauge  reactivation  via  a  second  substrate  and  do 
not  measure  OP  turnover  or  its  products  directly. 


Phosphonylated  esterases  may  either  spontaneously  reactivate 
or  reside  in  an  inactive  state,  which  is  the  aged  state  [8] .  However, 
for  wt  hCEl  no  aging  has  ever  been  measured.  As  shown  in 
Figure  IB  these  two  pathways  compete  as  first-order  reactions. 
Therefore,  the  true  rate  of  enzyme  dephosphonylation  (kp)  was 
determined  by  equation  2  [25] : 

*3  =  (4M*)(*«fa)/100  (2) 

where  kobs  and  Amax  were  defined  in  eq.  1.  The  denominator 
normalizes  for  Amax.  .411  data  were  analyzed  in  KaleidaGraph  4 
(Synergy  Software,  Reading,  PA).  Because  the  reactivation  step 
represented  by  k%  was  the  slow  step  in  the  OP  hydrolysis  reaction, 
ko,  effectively  represented  the  rate  of  OP  turnover. 

T o  determine  the  bimolecular  rate  of  inhibition  [kpj  with  racemic 
cyclosarin,  100  nM  of  purified  V146H/L363E  hCEl  or  wt 
enzyme  was  exposed  to  100-fold  excess  of  cyclosarin  and  the 
decrease  in  pNPB  hydrolysis  was  measured  over  time.  Data  were 
fit  to  equation  3  and  the  k\  was  calculated  by  dividing  the  rate  of 
phosphonylation  (kp)  by  the  cyclosarin  concentration: 

A=A0^t  (3) 

where  A  is  the  CE  activity  at  a  given  time  t,  and  A0  is  the  enzyme 
activity  prior  to  cyclosarin  inhibition. 

Further  inhibition  kinetics  were  determined  with  stereoisomers 
of  OP  model  compounds,  where  100  nM  of  V146H/L363E  hCEl 
was  incubated  at  room  temperature  with  increasing  concentrations 
of  inhibitor.  The  thiomethylated  OPs  examined  (Figure  1A)  have 
previously  been  shown  to  form  identical  adducts  to  bona  fide  nerve 
agents  [8, 1 3] .  These  compounds  were  generously  provided  by  Dr. 
John  Cashman  at  the  Human  BioMolecular  Research  Institute 
and  their  synthesis  is  described  by  Gilley  et  al.  [13].  Aliquots  of 
enzyme  in  the  presence  of  sarin,  soman,  and  cyclosarin  model 
compounds  were  removed  at  various  time  points  (up  to  1  hour) 
and  the  level  of  remaining  enzyme  activity  was  determined  by 
comparing  4-methylumbelliferyl  acetate  (4-MUA)  hydrolysis 
relative  to  an  uninhibited  sample.  4-Methylumbelliferone  fluores¬ 
cence  emission,  measured  at  450  nm  following  excitation  at 
350  nm,  was  monitored  at  37°C  on  a  Pherastar  microplate  reader 
(BMG  Labtech,  Cary,  NC)  and  the  data  were  fit  to  equation  4 
[26]: 

_At_  =  K;l  1  1 

Alnr  A'2  \OP\(\—a)  k 2 
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where  Km  was  the  nerve  agent  model  Michaelis-Menten  constant, 
ko  the  unimolecular  phosphonylation  rate  constant,  v  the 
remaining  percent  enzyme  activity,  and  [OP]  the  OP  concentra¬ 
tion.  The  term  a  was  defined  as  [S]  /(Am+[S]),  in  which  [S]  was  the 
substrate  concentration  and  Km  was  the  4-MUA  Michaelis- 
Menten  constant.  For  these  experiments  a  was  0.91.  All 
experiments  were  performed  in  triplicate  and  data  were  analyzed 
to  determine  fa,  where  ki  —  k^/K^  [16].  Again,  all  data  were 
analyzed  in  KaleidaGraph  4. 

Computational  Modeling  of  hCEl  Mutants 

To  simulate  the  hCEl-OP  crystal  structure  for  computational 
design,  a  new  residue  type  was  created  for  S221  covalendy 
inhibited  by  sarin,  soman,  or  cyclosarin.  Following  removal  of  all 
water  molecules,  saccharides,  ligands,  and  extra  monomers  from 
the  hCE  1 -cyclosarin  homotrimer  crystal  structure  (RSCB  PDB 
3K9B  [8]),  protein  hydrogen  atoms  were  added  into  one  monomer 
with  Molprobity  [27].  OP-inhibited  S221  hydrogen  atoms  were 
added  using  Accelrys  Discovery  Studio  Visualizer  2.5  (Accelrys, 
San  Diego,  CA).  Any  strain  or  steric  interactions  in  the  resulting 
structure  were  then  removed  by  performing  50  steps  of  conjugate 
gradient  energy  minimization. 

Residues  with  side  chains  pointing  towards  and  within  5  A  of 
the  OP,  VI 46,  or  L363  residues  (using  3K9B  numbering,  residues 
89,  90,  93,  96,  97,  100,  101,  145,  220,  222,  252,  254,  255,  304, 
318,  358,  359,  361,  364,  388,  425,  426,  and  468,)  were  allowed  to 
sample  alternative  conformations  during  the  design  phase  of  the 
simulation,  but  their  identities  were  not  modified.  At  positions  1 46 
and  363,  only  histidine  or  glutamic  acid  was  allowed,  respectively. 
A  standard  backbone-dependent  side  chain  rotamer  library  with 
expansions  by  one  standard  deviation  about  and  %2  was  used 
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